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ABSTRACT 


An experimental investigation was conducted to study convective heat transfer in 
Straight and curved rectangular channels of high aspect ratio that approximate plates 
of infinite extent. Experiments were performed at steady state in the turbulent flow 
regime with one wall held at a constant heat flux and the opposite wall essentially 
adiabatic. The effect of curvature induced secondary flow on heat transfer on the 
concave and convex walls was observed by comparing Nusselt numbers for four 
different configurations at several different Reynolds numbers. Significant heat transfer 
enhancement was observed on the concave wall. Correlations for Nusselt number as a 


function of Reynolds number were calculated for the cases studied. 
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total radiation resistance 
Taylor number 
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I. INTRODUCTION 


A. BACKGROUND 

Lord Rayleigh [Ref. 1] was seemingly the first to investigate the mechanics of 
flows with curved streamlines. His analvsis showed that inviscid curved streamline 
flows are stable providing that circulation increases with radial distance (for flow 
between concentric cylinders this condition 1s met when the outer cylinder spins more 


quickly than the inner cylinder). 


Later, Taylor [Ref. 2] studied viscous flow between concentric rotating cylinders 
and showed that when the inner cylinder is rotated while holding the outer cylinder 


Stationary, the onset of instability is predicted by the dimensionless parameter: 


Ud (2 
Ta (Taylor number)= — —-—_ 241.3, (eqn Tie 
v E 


i 


Where U is the tangential velocity of the inner cylinder, d is the gap between cylinders, 
r, is the radius of curvature of the inner cylinder, and v is kinematic viscosity. The 
instability is manifest in the formation of regularly spaced, counter-rotating vortices 
whose axes of rotation lie in the circumferential (main flow) direction (see Figure 1.1). 
Furthermore, Tavlor concluded that these vortices form in the laminar flow regime but 


do not induce turbulent main flow [Ref. 3]. 


The flow of a viscous fluid in a curved channel was studied by Dean [Ref. 4]. [le 
developed an analytical solution for flow in a curved channel approximated by 
concentric cylinders where the radius of the inner cylinder was much greater than the 
gap between cylinders. Instead of rotating the inner cylinder Dean imposed a pressure 
gradient along the channel wall and predicted that flow instability would occur when a 


critical value of a dimensionless parameter similiar to the Taylor number was exceeded: 


Ud vc 
De (Dean number)= — -—— = 36, (eqn: li 
V e 


i 


where d represents the channel half-width. 





Gortler [Ref. 5} went on to investigate viscous flow in a curved channel and 
concluded that the flow instability predicted by Dean would take a form analagous to 
the Taylor vortices of concentric cylinder flow (see Figure 1.2). Regularly spaced 
counter rotating vortices form along axes of rotation parallel to the main flow 
direction. Yih and Sangster [Ref. 6] concluded that the centrifugal force coupled with 
the velocity gradient near the concave wall caused instability due to a stratification in 


specific weight that led to vortex formation. 


The secondary flows predicted by Taylor and Gortler (Tavlor-Gortler vortices) 
have been confirmed experimentally using various flow visualization techniques, and 
hot wire anemometry [Ref. 7]. The presence of these vortices in the laminar flow 
regime led to speculation of the effect of this secondary flow on heat transfer over 
curved surfaces. Extensive research has been done which yields good experimental 
agreement With analytical solutions for laminar flow in curved channels of rectangular 
cross section [Refs. 8,9,10], and for turbulent flow in a curved rectangular channel of 
low aspect ratio (ratio of the spanwise dimension to the channel height) 
[Refs. 11,12,13]. Research in these areas indicates that there 1s enhancement of heat 
transfer in these flow regimes When Tavlor-Gortler vortices form near the concave wall. 
However, information relating to the fluid mechanics and heat transfer tn the transition 


and turbulent flow regimes for channels of high aspect ratio is very limited. 


B. APPLICATIONS 

Studying heat transfer over a curved surface is just one aspect of the complex 
fluid dynamics and heat transfer problems that exist in an operating gas turbine. The 
high tenyperatures characteristic of a combustion engine constrain gas turbine 
performance because of metalurgical hmutations of turbine blades which in turn 
restricts the efficiency of the turbine. Information gathered tn this” 1\peumem 
experimentation could lead to a better understanding of the problem and facilitate the 


design of gas turbines of higher efficieney and power to Weight ratio. 


CC. INTENT OF THIS STUDY 
The purpose of this investigation was to study the effect of curvature on heat 


transfer in a curved rectangular channel of large aspect ratio by gathering extensive 
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Figure 1.2 Taylor-Gortler Vortices in a Curved Channel. 
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data relating Nusselt number to Reynolds number for flow over walls of constant heat 
flux of straight, convex and concave geometry. This information will provide a firm 
data base for validating computer models of the complex flow and heat transfer in 


curved rectangular channels similar in geometry. 


Data for hydraulically fully developed flow in a straight rectangular channel 
heated either from above or below was taken to provide a basis for comparison with 
the flow over a heated convex or concave wall. Empirical correlations similiar to the 
familiar Dittus-Boelter correlation for Nusselt number as a function of Reynolds 
nuniber were calculated from the data so that the heat transfer enhancement could be 


described quantitatively and compared to previous data. 


A secondary goal was to determine the relative effect of buovancv on heat 
transfer in a curved channel. The effect of buoyancy on heat transfer in forced flow 
must be considered if the effects of free and forced convection are comparable because 
the combined heat transfer is then dependent on the direction of free convection in 


comparison to that of forced convection. 
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Il. EXPERIMENTAL APPARATUS AND PROCEDURE 


A detailed description of the materials and construction methods used to build 
the experimental channel and associated instrumentation 1s contained in G. Galvyo’s 
thesis [Ref. 14: pp. 25-43]. Some minor changes had to be made to the existing 


hardware and computer software and are described herein. 


A. DESCRIPTION OF THE APPARATUS 

The channel was built of plexiglass, having a channel width of 25.4 cm and a 
height of 0.635 cm, yielding an aspect ratio of 40:1. The channel was rigidly supported 
On a wooden frame and insulated over its entire length with three lavers of half inch 
foam insulation. The curved section subtends 180° of arc, and the convex wall (inner 
curved wall) had a radius of curvature of 29.7 cm (see Figure 2.1). Air was pulled 
through the channel by a compressor joined to the duct by two inch PVC pipe through 


a throttle valve. 


Four 30.48 cm x 25.4 cm wire wound rubber heaters laminated to copper plates 
Were mounted on the inside walls of the channel to accomplish heating of the air. Two 
heaters were mounted in the straight test section, one on the upper wall and one on the 
lower, at a sufficient distance from the channel entrance to insure that flow was 
hydrodynanucallv fully developed prior to flowing over the straight heated sections. 
The two remaining heaters were installed in the curved test section, one on the concave 
wall and one on the convex wall, with the leading edges of the heaters beginning at an 
angle of 101.5° from the curved section entrance; this allowed vortices to form prior to 
flowing over the curved heated sections. The heaters were constant power density type 
capable of 600 watts each, and could be energized individually or tn any combination 


of pairs by use of two power supplies. 


Thirty copper-constantan thermocouples were imbedded in each of the copper 
plates at the copper plate-rubber heater interface so that local temperatures could be 
measured (see Figure 2.2). Four thermocouples were mounted laterally across the 
channel entrance to measure inlet temperature (see Figure 2.1). Simuliarly. four 


thermocouples had been mounted at the exit of the channel; however, experimentation 
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Figure 2.2 Thermocouple Placement in the ITleated Plates. 
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30.48 cm 


showed that there was insufficient mixing of the air up to the exit which resulted in 
erratic outlet temperature measurements [Ref. 14: p. 67]. Instead, eighth inch 
honeycombing was installed at the exit of the channel to promote muxing, and the four 
exit thermocouples were placed in the PVC pipe near the end of the channelNeumay 
transition bell (see Figure 2.1). Two layers of half inch foam insulation were applied to 
the PVC to a distance of approximately fourty-eight centimeters past these 
thermocouples to minimize heat loss to the surroundings prior to measuring outlet 
temperature. One thermocouple was placed downstream of the flow measuring orifice 
and used to calculate the density of the air flowing through the orifice. Three groups 
of five thermocouples each were connected in parallel between the lavers of insulation 
surrounding the heated sections so that losses to the surroundings through the 


insulation could be estimated. 


Anv combination of eighty voltages from the thermocouples and power supplies 
could be read by the Hewlett Packard 3054A automatic data acquisition system used 
for this research. Coupled to the data acquisition system were a Hewlett Packard 9826 
computer and a printer which provided control of the data acquisition system and 


facilitated recording, reducing, and printing of data. 


A flow measuring orifice was installed in accordance with the ASME Power Test 
Code [Ref. 15], and water and mercury manometers were used to measure the pressure 
upstream of the orifice and differential pressure across the orifice so that mass flowrate 


could be caleulaved. 


B. EXPERIMENTAL PROCEDURE 

The procedures followed for testing were the same irrespective of the specific test 
configuration. The procedures are: 

1. Zero the manometers, shut the svstem throttle valve, then start the compressor. 
2. Use the throttle valve to adjust the differential pressure across themiiGm 
measuring orifice to that corresponding to the desired flowrate or Reynolds 
number. Determining the approximate differential pressure corresponding to a 
specific Reynolds number takes some experience and the use Of jams 


approximation: 


a) 


Res, » 
a) a bine | (AP), (eqn 2.1) 


3. Energize the heater and adjust the voltage for the approximate value of 


(Tout 
centigrade was used for the temperature differential, the objective being to 


- T. ) desired. In these experiments a value of six to twelve degrees 


maximize (T.,,-T.,) while linuting the highest local plate temperature to a 


out 
maximum of about 100°C. Maximizing this temperature differential lessens the 
uncertainty of the calculated Nusselt number. 

4. Experience showed that the channel reached steady state about four hours after 
first energizing the heater, and approximately one hour to steady out between 
consecutive data runs at different flowrates. When all of the following criteria 
Were met, steady state was considered to have been achieved: 


a. Change in the quantity (T 


g lee eGeeless than two percent over a ten 
minute period. 


out 


b. Change in the heated wall temperature of less than one percent over a ten 
nunute period. 


c. Change.in heated wall insulation temperatures of less than one percent over 
a ten nunute period. 


5. When the above criteria for steady state are met, run the data gathering and 
reducing program (see Appendix A) several times in succession to verify that 
Steady state has been reached. Three to six consecutive data collection runs at 
Steady state meeting the criteria of (4) above are required to determine an 
average Nusselt number versus Reynolds number relationship. Generally the 
run-to-run variation in Nusselt number and Reynolds number were much less at 
lower flowrates. 

6. When a sufficient number of consistent data runs have been recorded, throttle 
to adjust flowrate to the next desired Revnolds number and adjust the heater 


voltage to achieve the desired plate temperature. 


The atmospheric pressure, orifice differential pressure and pressure upstream of 
the orifice had to be read from a barometer and two manometers and entered into the 
data reduction program interactively; these parameters were not read by the automatic 


data acquisition system. 


During this study it was noted that small (one degree) ambient temperature 
perturbations had a noticeable effect on temperatures measured throughout the test 
system; therefore, close control of ambient temperature oscillations was required so 
that temperature changes in the test system could be mininuzed and consistent data 


measured from run to run. 


The initial configuration of the test channel was such that buoyancy induced flow 
opposed forced flow in the curved test section. After data was taken over the full 
range of flow that the compressor was capable of, the rig was inverted and data was 
taken with buoyancy induced flow aiding forced flow. The purpose of this aspect of 
the study was to determine whether natural convection contributed significantly to the 


total heat transfer in the curved test sections. 


rt 
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Hil. DATA REDUCTION 


Each heater plate provided a constant surface heat flux. The heater plates were 
energized by regulated D.C. power supplies through stable precision resistors of low 
resistance value (see Figure 3.1). The voltage drop across both the heater element and 
the precision resistor were read and recorded by the data acquisition system. Then the 


total power supplied by the heater was determined from the relation: 


Ven 
eho pr 


Q 
p ae 


feqn 3.1) 


where Ne was a Known value of resistance included as a constant in the data reduction 
program. The insulated, unheated walls were treated as adiabatic because the heat loss 
through them was negligible compared to the total heat input from the heater plate. 
Therefore, the bulk temperature varied linearly from the leading to the trailing edge of 


the heated plate. 


Thermocouple voltages were converted to centigrade temperature using individual 
thermocouple calibration curves. All temperatures were recorded then combined to 


calculate specific average temperatures desired. 


The actual heat loss to the surroundings was difficult to account for. Typically 
the estimated total heat loss by conduction through the insulation and by radiation 
was only about ten percent of the actual total heat loss. The major difficulty in 
making an accurate prediction of the heat loss was due to not being able to determine 
the area of the channel over which heat was lost to the surroundings. For the purpose 
of estimating the heat loss, the area of the heated plate was used as the effective area 
of loss: however, the plexiglass section of the channel adjacent to the heater plates acts 
as extended heat transfer surfaces, and the other copper heater surfaces in the channel, 
which have a low resistance to thermal conduction, act like heat sinks. Therefore, 
instead of assuming that the energy transferred to the air was equal to the energy input 
from the plate minus the estimated energy lost to the surroundings over the area of the 


plate, the energy actually transferred to the working fluid was calculated as follows: 


qd 
7 Heater 


Vp “@—Precision Resistor 


Power Supply 


Figure 3.1 Power Measuring Circuit. 


Qair = Me, AT (eqn 3.2) 
umcnce Wee LoL. 
The corresponding plate average heat transfer coefficient and Nusselt number based on 


the hydraulic diameter of the duct were then calculated as follows: 





h = 22 (eqn 3.3) 
Ant (Ty - Ty) 

_ hD 

“Ui ai (eqn 3.4) 


Where D,,, the hydraulic diameter of the channel, is a geometric constant, 


4 X Cross Sectional Area of the Duct 


D eqn 3.5 
hd Wetted Perimeter of the Duct 5 


Data on atmospheric pressure, orifice differential pressure, and orifice upstream 
pressure was used with the flow measuring orifice characteristics and air densitv based 
On air temperature at the orifice to calculate mass flowrate, m. Reynolds number 


based on hydraulic diameter was then found by the relation: 





Re, = —23 (eqn 3.6) 


where [U... was the viscosity of the working fluid at the average bulk temperature 
measured in the test channel. Similiarly, Reynolds numbers based on the pipe diameter 
and channel height were also determined. For the concave and convex sections only, 


the dimensionless parameter, Dean number, was calculated from: 


a m D /2 D /2 
Aft ie 


ce" air 1 





(eqn 3.7) 
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IV. EXPERIMENTAL RESULTS 


Several sets of data were taken over the approximate range of Reynolds numbers 
from 6000 to 20,000. A computer program (see Appendix A) was used to reduce the 
raw data and calculate the quantities of interest, particularly plate average Nusselt 
numbers for corresponding Reynolds numbers, both based on the hydraulic diameter of 
the channel. Three to six separate data sets at each Revnolds number were taken, then 


averaged to give the value of the quantities of interest. 


After extensive data had been taken for the channel configured such that free 
convection opposed the direction of forced convection, the channel was inverted and 
several sets of data were taken for the convex and concave plates with the directions of 
free and forced convection coincident so that the effect of buovant flow could be 


deternuned. 


The significant parameters and corresponding uncertainty for the different 
configurations follows in Tables | through 4. Graphic comparisons of plate averaged 
Nusselt number versus Revnolds number are shown in Figures 4.1 through 4.7. These 
graphs also include calculated uncertainty ranges for the data, the best fit correlations 
for each configuration, data for comparison from other sources, and some of the 


correlations that have been documented in research similiar to this experimentation. 


A. HEATING THE STRAIGHT UPPER PLATE 

A total of sixty-nine data pairs were taken over an approximate range of 
Revnolds numbers of 6800 to 20,000 (see Table 1). One means used to analyze the 
data was the statistics software, MINITAB [Ref. 16]. It was observed that the data at 
Revnolds numbers below 7000 exhibited a relativelv large deviation from the regression 
in comparison to the data for Reynolds numbers above 7000. This result agrees with 
previous research done by J. Wilson [Ref. 17: p. 50] which revealed that transition flow 
in the channel persisted up to a Reynolds number of approximately 7100; therefore, 
data at Revnolds numbers below 7000 was not included in the regression calculated for 


the upper straight plate. 


Using a least squares regression, a correlation for the heat transfer from the 


upper plate was determined to be: 
Ma, —wro2dRe, (eqn 4.1) 


where Reynolds and Nusselt numbers are based on the hydraulic diameter. The 
Statistical correlation of the regression was 0.998. The experimental uncertainty of the 


data was approximately £5% for Nusselt number and £6% for Reynolds number. 


For comparison, consider the well Known Dittus-Boelter correlation based on 


hydraulic diameter and a Prandtl number of 0.71: 
Nu, = 0.02Re, 7° (eqn 4.2) 


At first glance this small difference between the correlations seems insignificant; 
however, a look at Figure 4.1 shows that the Dittus-Boelter correlation overestimates 
the actual Nusselt number by about thirty to thirty-seven percent over its range of 
applicability (Re 2 10,000). 


Also shown for comparison on Figure 4.1 are interpolated values from the Kays 
and Leung [Ref. 18] analytical solution for infinite parallel plates. These values too lie 
outside the uncertainty limits of the actual data, but only overestimate Nusselt number 


by ten to fifteen percent when compared to the finite flat plates of this research. 


Data previously taken by G. Galvo [Ref. 14: p. 58] 1s included in this plot. His 
data agrees with that of this research to within the bounds of uncertainty. It must be 
noted that the data taken by Galyo at Reynolds numbers above about 17,000 1s less 
certain than the other data due to the difficulty he encountered in recording an 


accurate channel outlet temperature at high flowrates [Ref. 14: p. 67]. 


B. HEATING THE STRAIGHT LOWER PLATE 

A total of sixtv-four data pairs were recorded over an approximate range of 
Reynolds numbers of 6000 to 19000 (see Table 2). As in the previous case for the 
upper plate, the data recorded for Reynolds numbers less than about 7000 were not 
included in the least squares regression because of their relatively high deviation from 


the regression. The best fit regression correlation was: 


pe) 


NG | SO UlenRe ea (eqn 4.3) 


The statistical correlation of the regression was 0.998, and the experimental uncertainty 


of the data was essentially the same as that for the upper flat plate. 


Figure 4.2 shows that over the full range of testing there is essentially no 
difference in the results for the two flat plate cases. Data previously taken by G. 
Galyo [Ref. 14: p. 59] is shown as well for comparison, and is well within the 
experimental uncertainty of the data recorded in this study except at the highest 
Revnolds number where, as mentioned previously, G. Galyo’s data was highly 


uncertain. 


C. HEATING THE OUTER CURVED (CONCAVE) PLATE 
Seventy-two data pairs were taken over an approximate range of Reynolds 
numbers of 6000 to 20,000 (see Table 3). A least squares regression of the data 


resulted in the relation: 
Nii, = O05Ree 2: (eqn 4.4) 


The statistical correlation of the regression was 0.999, and the experimental uncertainty 


was the same as that for the two previous cases. 


Figure 4.3 revealed that the previous data reported by G. Galvo [Ref. 14: p. 60] 
coincides well with the data recorded in this study except, again, at the highest 
Revnolds number where G. Galvo had difficulty determining duct outlet temperature 
accurately. Also shown for comparison are extrapolated correlations from Kreith 
[Ref. 12: p. 1253]. and Brinich and Graham Ret. apace 


[In contrast to the results for the flat plates there was no significant deviation 
from the regression even at the lowest Revnolds number tested for the concave plate. 
This supports previous research bv J. Wilson [Ref. 17: p. 50] which indicated that the 
limiting Reynolds number for transition flow in curved sections was less than that for 


flat plates and approximately equal to 3300. 


The data for the flat plates was used as a baseline for comparison with the 
concave plate, and heat transfer enhancement for the concave section was reflected by 
an increase in Nusselt number of about twenty-one percent at a Reynolds number of 


10,000 to fourteen percent at a Reynolds number of 20,000 (see Figure 4.4). 


There is a tendency to extrapolate this data to where the correlations for the 
concave plate and the flat plates converge in order to estimate a Reynolds number 
where the effect of curvature is no longer significant. For these correlations 
convergence occurs at a Reynolds number of about 205,000 and a Nusselt number of 
approximately 232; however, research by Brinich and Graham indicated that the 
enhancement in heat transfer in a concave section was significant even at Reynolds 
numbers as high as 250,000. 


D. HEATING THE INNER CURVED (CONVEX) PLATE 
A total of seventy-three data pairs were recorded over a range of Reynolds 
numbers of 6000 to 20,000 (see Table 4). A least squares fit of the data gave the 


relation: 
Rie= OOlRe, 7° (eqn 4.5) 


at a statistical correlation of 0.997, and an experimental uncertainty the same as for the 


previous cases. 


Figure 4.5 revealed that, as was the case for the concave plate, there was no 
significant deviation from the regression even at the lowest Revnolds number tested. A 
comparison of this data with that reported by G. Galvo [Ref. 14: p. 61] showed that 
there is agreement to within the bounds of experimental uncertainty except at the 
highest Revnolds number. For comparison, data reported by Kreith [Ref. 12: p. 1253], 
and Brinich and Graham [Ref. 11: p. 36] for a convex surface were included in the 


graph. 


Using the flat plate data as a reference of comparison, the reduction in heat 
transfer on the convex plate was manifest in a Nusselt number consistently six percent 
lower than that of the flat plates over the entire range of flow rates. This did not 


represent a large difference between the convex and flat plate cases: however, six 


2» 


percent was at least one percent more than the conservatively calculated experimental 
uncertainty of the data. Therefore, this data implied that a decrease in heat transfer on 


a convex surface should be expected, though probably small. 


A comparison of the convex and concave surfaces (see Figure 4.6) pointed out 
that the concave surface heat transfer is approximately thirty percent greater than that 
for the convex surface at a Reynolds number of $000, and about twenty percent greater 
at a Revnolds number of 20,000. 


E: EFFECT OF BUOYANCY 
Tables 3 and 4 include data taken when the test channel was oriented such that 
free convection aided forced convection. Graphical comparison for this case is shown 


in Figure 4.7. 


Data taken for this configuration was within the bounds of experimental 
uncertainty, and showed no pattern in differing from the case where the direction of 
free convection opposed that of forced convection. These results indicated that there 


was no significant buovant effect over the full range of flow rates studied. 
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* Indicates measurements that were taken With the test channel inverted 
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Figure 4.2 Straight Lower Plate Results. 
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Figure 4.3 Concave Plate Results. 
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Figure 4.4 Comparison of Straight Upper and Concave Wall Results. 
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Figure 4.5 Convex Plate Results. 
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V. DISCUSSION 


A. FLAT PLATES 
The experimental results for the upper and lower straight sections of the channel 


were essentially the same, and are not considered separately in this discussion. 


The Dhittus-Boelter correlation for Nusselt number based on the hydraulic 
diameter of the test channel and a Prandtl number of 0.71 yielded Nusselt numbers 
thirty to thirty-seven percent higher than those recorded in this experimentation. The 
Dittus-Boelter correlation is an empirical approximation that was developed to 
estimate local Nusselt numbers for circular tubes. The correlation itself is only an 
approximation, and when modified for non-circular ducts where convection coefficients 
vary about the periphery, the approximation becomes less certain. Here the correlation 
Was used to estimate the average Nusselt number over the length of the heated section. 
To a first approximation the Dittus-Boelter correlation can be used to estimate the 
average Nusselt number because the ratio of the plate length to the hydraulic radius 
Was greater than ten: however, this assumption adds to the uncertainty of the resulting 


convection coefficient. 


The solution for parrallel plates of infinite extent developed by Kays and Leung 
gave Nusselt numbers closer to the experimental results, but was still high by ten to 
fifteen percent, though ten to fifteen percent accuracy when estimating heat transfer 
coefficients in the turbulent flow regime is not unusual. Even though the straight 
section of the test channel approximated flat parallel plates of infinite extent, some 
sidewall influence on flow was expected and mav account for a part of the disparity of 
the comparison. The graphical comparison shown in Figure 4.1 was developed by 
linear interpolation of the data from [Ref. 18: p. 552], though there was no evidence to 
support the assumption that the results of the Kays and Leung solution are a linear 
function of either Prandtl or Revnolds number; therefore, additional inaccuracy may 


have been introduced into the comparison. 


B. CURVED PLATES 

Experimental and analytical investigation of flow and heat transfer on concave 
and convex surfaces in a rectangular duct of high aspect ratio has thus far received 
little attention; consequently, there is scant previous work available for comparison. 
Some results are available for channels of modest aspect ratios, but these do not 
generally compare well due to the significant sidewall effect introduced at low aspect 


ratios. 


Studies of the flow and heat transfer characteristics in high aspect ratio channels 
has been done at the Naval Postgraduate School [Refs. 9,14,17,19,20,21,22]. Results of 
this previous research were comparable to the results of this experimentation and were 
highly depended upon to validate these results. Only Galvo’s work 1s shown for 
comparison in this report, but his results agreed within the limits of experimental 
uncertainty with those of Wilson [{Ref. [7], and Daughety [Ref. 22]. The major 
Improvements of this experimentation over that previously done were the decreased 
experimental uncertainty, especially at the highest Reynolds number, and _ the 
preponderance of data that permitted empirical correlations to be determined to a very 


high level of confidence. 


fet) Kreith | Ret. 12: p. 1253) investigated the influence of curvature on heat 
transfer in a rectangular channel of aspect ratio equal to three. For a Prandtl] number 
of 0.71, graphical comparison with the results recorded in this experimentation showed 
that Kreith’s results were approximately thirty percent lower (see Figures 4.3 and 4.5); 
however, Kreith’s test apparatus not only had a much lower aspect ratio, it was also 
heated by a heating element that was centered between the upper and lower walls of 


the channel instead of heating the walls themselves. 


Research done in 1977 with a channel of aspect ratio six by Brinich and Graham 
[Ref. 11: pp. 35-36] was shown for comparison in Figures 4.3 and 4.5. The range of 
Reynolds numbers over which they performed their experiments was approximately 
$4,000 to 253,000. For the concave surface, Brinich’s and Graham's data extrapolated 
to a Reynolds number of 15,000 vielded a Nusselt number about fiftv percent higher 
than that recorded in this experiment. Likewise, an extrapolation to a Reynolds 
number of 15,000 for the convex plate showed that there was about a thirty percent 


difference in the results of the two experiments. These differences are rather large, but, 
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again, the two studies differed greatly in that the aspect ratios and the range of 


Reynolds numbers investigated were much different. 


Similiar research was done in 1978 by Uchida and Koizumi [Ref. 13: pp. 
1714-1716], with channels of aspect ratio up to four; however, the full report has not 
been translated into English. From figures, graphs, and equations included in their 
report it appears that their test channel was heated along its entire periphery and that 
the published results were in the form of an average Nusselt number for the 


combination of a convex and concave wall. 


It is unlikely that an average heat transfer coefficient or Nusselt number can be 
deternuned when both the convex and concave walls are heated by simple 
superposition of the individual correlations for the convex and concave walls because 
of the complex nature of the flow in a curved channel. Therefore, there was no 
manner in which the results detailed in Ref. 13 could justifiably be compared to the 


results of this experimentation. 


Even though the test equipment was such that a pair of heaters could be 
energized to supply a constant heat flux on opposite walls of the channel, no 
conclusive experiments could be performed to determine an average Nusselt number 


for this configuration. 


Referring to (eqn. 3.3) and denoting hy as the average heat transfer coefficient for 
single plate heating, and h, as the average heat transfer coefficient when heating a pair 


of opposite walls. then: 


(a te 


h ef Wo Wi : Ton 


1, (T 


at 


(eqn 5.1) 


— 


wh - Tot? 


Where T.,, represents the average temperature when one wall (inner or outer) is heated, 
and the heat flux from both the inner and outer plates was assumed equal when 
heating a pair of plates. The average wall temperature term in the numerator of (eqn. 
5-1) cate De Tecast.as- 


i eee 
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where OT represents the difference between the average temperature when heating a 
pair of plates and the single plate average temperature. Combining (eqn. 5.1) and (eqn. 
ee) yields: 


rn : Tank + oT) (eqn 5.3) 
We i: Totk) 
1+ 8T(T yy - Ty) 
l 
1 + oT/AT 


l 


Pus 


(Eqn. 5.3) reveals that only when the temperatures of both plates are the same, such 
that OT = 0, when heating a pair of plates, can a true average heat transfer coefficient 
for both walls be obtained. Strictly speaking the temperature at all corresponding 
points opposite each other on the plates should be equal to vield an average heat 
transfer coefficient for the combination of both walls. The heaters used in this 
experimentation were constant heat flux heaters, and teniperatures at corresponding 
points on opposite walls could not be matched; therefore, conclusive data for the case 


of heating a pair of plates could not be gathered using the existing equipment. 


Though local heat transfer coefficients and Nusselt numbers were calculated at 
the leading edge, at mid-plate, and at the trailing edge of the heated plates, no 
conclusive pattern was evident in the data. The same was true for local temperatures. 
More conclusive results could not be obtained because there were too few 
thermocouples spaced too far apart; that 1s, the spanwise thermocouple spacing was 


greater than the diameter of the Tavlor-Gortler vortices that formed. 


VI. CONCLUSION 


The results of this experimentation, detailed in Chapter IV, indicated that 
secondary flow in a curved rectangular channel in the form of Taylor-Gortler vortices 
influenced the heat transfer on the concave and convex surfaces over a range of flow 


corresponding to Reynolds numbers from about 6000 to 20000. 


Correlations of Nusselt number as a function of Reynolds number for the two 
flat plates, the convex plate, and the concave plate were calculated to a high level of 
certainty and showed that heat transfer was enhanced by approximately twenty percent 
on the concave surface as compared to the flat surfaces. Similiarly, heat transfer over 
the convex surface was reduced by approximately six percent in comparison to the flat 
surfaces. Furthermore, the enhancement or reduction in heat transfer on the curved 
plates was noted to be greatest at the lowest Revnolds number tested, and tapered off 


with increased flow rate. 


Comparisons made of the results obtained in this experimentation to that of 
other pertinent research supported the conclusion that heat transfer is enhanced on a 
concave surface. Experiments performed for the purpose of determining the relative 
significance of buovancy induced flow indicated that even at the lowest Reynolds 
number of 6000, buoyant free convection was negligible compared to forced 


convection. 
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VII. RECOMMENDATIONS 


EQUIPMENT 


Small oscillations in ambient temperature (+1°) noticeably affected 
temperatures in the channel and contributed to the differences in the results in 
successive data taking runs. Maintaining close control over ambient 
temperature would make experimental results more reproducible and decrease 
the time necessary to gather data. 

Existing equipment allowed for testing onlv up to a Reynolds number of about 
20000, though the corresponding Revnolds numbers for flow through a gas 
turbine are much higher. A larger test channel and higher capacity air 
compressor would be required for experimentation at higher flow rates. 

A great amount of time was spent waiting approximately four hours for the test 
channel to reach steady state after first energizing the heater. Installing guard 
heaters would significantly reduce the warm up time to. steady state. 
Additionally, guard heaters would prevent heat loss through the insulation 
surrounding the heaters. 

Pressure sensors that could be read directly by the automatic data acquisition 
system several times during the course of a data run would increase the 
accuracy of flow rate measurements and permit successive data runs to be taken 
without the presence of an operator. 

More thermocouples should be installed closer together to improve the accuracy 
of temperature measurements and to facilitate studving the pattern of local 
variations in temperature and heat transfer in the spanwise and flow-wise 
directions. 

The single most significant contributor to experimental uncertainty in the 
calculated average Nusselt number and hydraulic Reynolds number was the 
uncertainty in the channel cross section. The uncertaintv in experimental 
results could be greatly reduced if a channel were constructed paying special 


attention to the cross section dimensions. 
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B. 


FURTHER RESEARCH 


Experimentation similiar to this but at higher flow rates should be done so that 
a large data base could be established and heat transfer correlations could be 
calculated for a wide range of application. 

Experiments at lower flow rates where buoyancy 1s significant should be done 
so that a suitable length parameter can be deduced for calculating the Grashof 
number for curved channel flow. 

Different curvatures should be studied to determine the effect of the degree of 
curvature on heat transfer in a curved channel. 

Thermostatic control should be included in the experimental equipment that 
could maintain equal plate average temperatures when heating a pair of heaters 
on opposite walls so that experimental results for channel average heat transfer 


coefficients could be gathered. 
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APPENDIX A 
DATA REDUCTION PROGRANI 


ee lee NAME: MAIN 
EaoleReEVISION DATE: 19 OCT 1986 
BD EIEDEEY LisedoR. HANK II] 


! PROGRAM FOR GATHERING AND REDUCING DATA WHEN HEATING ONE COPPER PLATE, 
! WITH S MINUTES BETWEEN RUNS. PROGRAM IS MODIFIED TO ACCOUNT FOR 

f BROKEN THERHMOCOUPLE CHANNELS S, 16, 21, AND 26 WHICH EFFECT OUTER 

! CURVED HEATER, THERMOCOUPLE CHANNEL 47 WHICH EFFECTS INNER CURVED 

! HEATER, AND THERMOCOUPLE CHANNEL 11 WHICH EFFECTS OUTER STRAIGHT HEATER 


COM /Co/ D(7) ,Aa(76,2) 
Meme 173) ,1(79) 


Pec OR nee LON TeACTORS FOR EMF 10 DEGREES CELSIUS FOR CU-CO I1/C. 


DATA 0.10086091 ,25727 .94369, -767345.8295, 78025595.81 

DATA -9247486589,6.97688E+11,-2.66192E+13,3.94078E+14 

READ DC #) 

PRINTER IS 701 

BEEP 

INPUT “ENTER MONTH, DAY, AND TIME (MM: DD:HH:NN:SS)",Times 

OUTPUT 709;"TD"; Times 

PEEL 

INPUT "ENTER INPUT MODE (0=3054A,1=FILE)*.Im 

IF Im=1 THEN 

BEEP 

INPUT “ENTER NAME OF EXISTING DATA FILE" ,Oldfiles 

PRINT USING “20X,""These Results Are From Data File: “",10A";Oldfirtle$ 
PRINT 

PRINT 

Boe 

BEEP 

INPUT “ENTER NAME OF NEW DATA FILE?" ,Newf i les 

PRINT USING "20X,""Data For This Run Is Stored In Data File: “".10A" 3Neuf 


PRINT 
PRINT 
END IF 


f ENTER THE CHANNEL CONFIGURATION. 


aa USING "9X,""THE FOLLOWING CHANNEL CONFIGURATION HAS SELECTED: """ 
RINT 

BEEP 

INPUT “SELECT HEATER TYPE (O=STRAIGHT, 1=CURVED)",Itype 

IF Itype=0 THEN PRINT WSING "12X,""“Heating STRAIGHT Test Section.""" 
IF Itype=1 THEN PRINT USING "12X,""Heating CURVED Test Section.""" 
BEEP 

INPUT “SELECT HEATER POSITION (0=OUTER, 1=INNER)",Ipos 

IEtpes=0 THEN PRINT USING “12X%.° "Heating OUIER Plate.""" 

aaa MHENSPRINTMUSTHG= "12x. Heating INNER Plate. << 
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2040 


1 ENTERS FILES WITH T/C CALIBRATION COEFFICIENTS BASED ON HEATER TYPE. 


IF Itype=0 THEN ASSIGN ®Filec ie) secs Ge lels y 
IF Itype=1 THEN ASSIGN @Filec 10 “CURV_COEF™ 
FOR I=0 TO 76 

ENTER aE ileciAa(1,0),AatI,1) Aall.2) 

NEXT I! 

ASSIGN 8Filec 1 Oe s 


1 INSERT DATA STORAGE DISK. 


BEEP 

INPUT “CHANGE DISK AND HIT ENTER , Ok 
IF Im=t THEN ASSIGN @File TO Oldfiled 
IF Im=0 THEN 

CREATE BDAT Newfirlesd.40 

ASSIGN @File TO Newfiled 

END IF 

ENTER 709; Timed 

CLEAR 709 

PRINT USING "90X, "Month, Day. and Time. °* 14A"; Times 
PRINT 


1 FOR DATA RUN READ AND STORE ALL RAW EMF VALUES. 


IF Im= THEN 

OUTPUT 709;"AR AFOO AL 79° 
AUTPUN reeset lan lteter 
FOR I=0 TO 79 

QguTPUT 709; "AS* 

ENTER 722;Emf (1) 

NEXT I 

OUTPUT @FilesEmt(*) 

Ee 

ENTER @Firle:Emf (*) 

END If 

guTPUT 709; "TD" 


1 CONVERT AND CALIBRATE ALL T/C READINGS TO DEGREES CELSIUS. 


PRINT USING eOx .° THE FOLLOWING DATA WERE RECORDED: °*" 
PRINT 

FOR [=0 TO 79 

ie 1=38.0R 1-32 OR I+77 OR 1278 THEN 

IF (Ipos20 AND [=38) OR (Ipos=1 AND 1=77) THEN 


Tadd=0 

IF Ipos#! THEN ladd=39 

PRINT USING 1 OX. "Heater Voltage (Vh) « )S 20 acne cy" sEmf (38+la 
PRINT USING "12X.""Resistor Voltage (Vr) = > ee cy" *"sEmfi(33+La 


IF Ipos=¥ THEN PRINT USING 41% “UP pec io hou Resistor (Rer?) * 9.01729 (ng 
IF Ipos=! THEN PRINT USING H1ay  “"Precis tom Rasistor «elms 2.00839 (Ohm 
PRINT 

END IF 

ELSE 


1 IF T/C QUT OF RANGE PRINT A QUE. 
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MO p anro—-OW 


Zo 20 
ya 


2530 


Aemt =Emft (1) 

IF Aemf<7.09E-4 OR Aemf>3.630E-3 THEN 
fel) =399.99 

Bese 

CALL Tyvsv(Emf(1I),Tt) 
IF I=79 THEN 

TeLy=Tt 

GOTO 2210 

END IF 

NCL) =FNieorr (lt ab) 
END IF 

END IF 

IER I 


! PRINT ALL T/C READINGS. 


USING t2x, CUTER PEATE: 
USING “14X,""T/C Number: 1 2 e: 4 5 


USING Seo UeOD ek). (CO Co C2 CSc C4), T 
USING 10 11 

USING 
US ING 
USING 
USING 
USING 


PRINT 
Boil 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


oiKe temp (6): 

"14X,""T/C Number: 7 8 3 

Caner Dep ex) 22116), 1079, 008). 609). 1 C10), 
13 iS 16 17 

eh E@epmppeeX) “:1(12) TS) , 114), 7015) , TC 
19 20 21 22 23 

ECON = 111), 1019), 1 (202) ), 1( 


Sake emp (Co 
"14%." "T/C Number: 14 


14k. temp (Car 


m7?) 


PRINT 
PRINT 


"“14X,""T/C Number: 
eva.  lemp (Co). 


ee), 1(23) 


2340 
2 liad 
2350 


26 26 27 28 29 
aS (ane mme ie 624) ,7(25),1(26), 127), 1 


USING 
USING 


PRINT 
PRINT 


"4X,  '7C Number: 
7 \4X Teme ve). 


Za),'(29) 


PRINT 
' CALCULATE AVERAGE TEMPERATURE OF OUTER HEATED WALL. 


Sum=0. 

IF Itype=0 THEN 
FOR I=0 TO 29 

IF [=11 THEN 2440 
Sum=Sum+t+T (1) 

NEXT I 

Two=Sum/29 


COR emo Oke BROKEN T/C CHANNELS 5S, 16, 21, AND 26 IN OUTER CURVED HE 


pioe 

FOR 1=0' TO 29 

IF [=S5 OR I=16 OR J=21 OR I[=26 THEN 2530 
Sum=Sum+T (1) 

NEXT I 

Two=Sum/26 

END IF 


PRINT USING *12X,""OUTER INSULATION: "*" 


= 


2570 PRINT USING °14X,°“T/C Number: 31 (liaise!) 32 (Tinso2d) 
2580 
2530 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
ae Sos (0. 
2700 ! 
2710! 
2720 
2130 
2740 
2750 
2/760 
2770 
2780 


PRINT USING “14X,“"Temp (C): 
Tinsol=T(39) 

Tinsoz=(C31) 

linsoc=10s2) 

PRINT 

PRINT USING “12X,”"°“ORIFICE TEMP (loert)) =) 330 0DDe 
PRINT 

PRINT USING “12X."“ INLET TEMPERATURE: *** 
PRINT USING °14X,""T/C Number: a5 
PRINT USING “14X,"*"*Temp (C): 
PRINT 


CALCULATE AVERAGE INLET TEMPERATURE. 


(C) "iia? 


36 37 Son 


Sum=Q. 

FOR I=34 TO 37 

Sum=Sumt+T (1) 

NEXT I 

Tin=Sum/4 

PRINT USING “24s "] NNER REATE: = 
PRINT USING “°14X,"“T/C Number: 4} 

46" a «@ 


2790 PRINT 
44) ,T(45) 
2800 PRINT 


2610 PRING 
S004 151) 
PRINT 


esa0. PRINT 
Sots.) 
2840 PRINT 
64"" u 

2850 PRINT 
62) ,1T(63) 
2860 PRINT 
7 fg te 

2870 PRINT 
638),T(69) 
2880 PRINT 
2390! 
Sa ! CALCULATE AVERAGE TEMPERATURE OF INNER HEATED WALL. 
2320 
2930 
294i} 
2oo0 
2960 
2500 
2980! 
23930 | 
3000! 
300s BoE 
3020 FOR [+40 TO 69 
3030 IF I[=47 THEN 3050 
3040 Sum=Sum+T(T) 

3050 NEXT I 


42 43 ad 
USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 


USING 


“lax,” Temp (Cy: 


"14X,"°"T/C Number: 47 48 49 50 


"14X,""Temp (C); 


"14X,""T/C Number: 51S) 54 SS 56 


“t4X,""Temp (CC); 


~14X,5 17C Number: 53 60 61 62 


"“14X,"*"*Temp (C): 


P1aXk.-17C Number: 65 66 67 68 


"Tak, Temps oC) | 


Sum=0Q. 

IF [type=0 THEN 
FOR [=40 TO 69 
Sum=Sum+T( 1) 
jew 

Twi =Sum/30 


CORRECTION FOR BROKEN T/C CHANNEL 47 IN INNER CURVED HEATER. 


“" .4(3D.0D.2X) 31134), [Cea 


33 (Tinsoaae 


T(37) 


45 


"",6(30.0D.2X)":7T(40) .T(41), 1042), | Cae 


51 


"™ 6(3D.DD,2X) "31(46).1(47),1(48), T( 49), 


S7 


"",6(3D.DD.2X)";7T(52),7(S3) , 1(54), hee 


63 


“* 6(30.D0D,2X)":7(58) .T(S9), T(60)). fcc 


69 


“* 6(3D.DD.2X)":1(64).1(65) ,7065) , 1 Cave 


3060 Twi2Sum/29 
3070 END IF 
3080 PRINT USING “12X."“INNER INSULATION: “"* 
SUOUee ERIN, USING@ET4X,"297C Number : viet linsit) 72 Ulinsi2) 73 Glens i 3) 
3100 PRINT USING “14X,"““Temp (C): WACO eOU oN) oll 70). br ti 1 Cie) 
eno Tinsihent (70) 
ei2z0) fTinsi2=s T(/1} 
es Tinsa3g=1(72) 
3140 PRINT 
eioO PRINI USING *12X,""QUILET TEMPERATURE:**” 
3160 PRINT USING "14X,""T/C Number: 74 75 76 aed 
e706 PRINT USING "4X, °"Temp (C): Peas elueen) + 1073), 0074) 1 (75). 1076) 
3180 PRINT USING "@,4" 
SSIOR. 
3200 ! CALCULATE AVERAGE OUTLET TEMPERATURE. 
e210)! 
3220 §=Sum=0. 
e2o0 FOR [=73 TO 76 
3240 Sum=Sum+T(T) 
3250 NEXT I 
3260 Tout=Sum/4 
e270! 
ae ' CALCULATE DIFFERENCE BETWEEN CHANNEL INLET AND OUTLET TEMPERATURES. 
S300 Idiff=Tout-Tin 
Bo10! 
3320 ! CALCULATE CHANNEL BULK TEMPERATURE. 
3330! 
3340 Tbhik=(TintTout)/2 
3350! 
on !' PRINT CALCULATED TEMPERATURES. 
3380 PRINT USING "9X,""THE FOLLOWING TEMPERATURES WERE CALCULATED:""* 
3390 PRINT 
3400 PRINT USING “12X,"*"Average Outer Wall Temperature (Two) = re. 
00 aia He)” “= Two 
3410 PRINT USING “12X,"“"Averagqe Inner Wall Temperature (Twi) Sle i 
epee (C)°""? Ti 
3420 PRINT USING “{2X,"“"Average Outlet Temperature (Tout) hedaach reo] B 
fee 6 6(C) Ys Tout 
3430 PRINT USING “12X,"““Average Inlet Temperature (Tin) = es D 
pi,” “ Me) ae Teh 3 
3440° PRINT’ USING “12X,""Channel Inlet and Outlet Temp Diff (Tdiff) ee |) 
epee «(C)""" sTdit ft 
3450 PRINT USING “12X,."“"°Averaqge Bulk Temperature (Tblk) ae eG) 
On. «€6CC Tbk 
3460! 
oe ' CALCULATE DELTA T (Tdel) FOR USE IN CALCULATING H LATER. 
‘> 
3490 Tdelo=Two-Tblk 
Bevu) Ideli=lwi-Tblk 
3510 Ir Ipos=0 THEN Tdel=Tdelo 
Sis IF Ipos=1 THEN Tdel=Tdeli 
3530 PRINT USING "12X,""Mean Temperature Difference (Tdel) ae Oe 
fee (C)"""* Tdel 
3540! 
aoe ! CALCULATE DIFFERENCE BETWEEN OUTER AND INNER WALL. 
3570) = Twdiff=ABS(Two-Tu1) 


53 


3960. 


PRINT USING "12X."“Quter and Inner Wall Temp Difference (lTudift)) jeueaeeeee 
CC Ot ereaincgl a tat 


! CALCULATE LOCAL CHANNEL TEMPERATURES AT INLET, MIDDLE AND OUTLET. 


Téern= linet S* lanehte eo 

Temid=Tint(4.5*Tdiff/12) 

feoutsslient Cll. o* lcdwety ico 

loos ee "12K," "Average Local Channel Temp at Inlet (Tein) = UE 
(Ce Hein 

PRINT eee "12X,"“Average Local Channel Temp at Middle (Tcmid) ee 20) 6 
(0) 8" enua 

PRINT USING “12X,“"“Averaae Local Channel Temp at Outlet (Tcout) eT WOUr 
CC)" ==" heout 


! CALCULATE AVERAGE LOCAL HEATED WALL TEMPERATURE AT INLET. 


Sum=0. 
Tadd=0 
IF Ipos=!I THEN Iadd=40 


! CORRECTION FOR BROKEN T/C CHANNELS 5S, AND 47 CURVED HEATER ONLY. 


IF Itype=1 THEN 

FOR I=0 TO 6 

IF (Ipos=0 AND I=S) OR (Ipos=! AND I=6) THEN 3810 
Sum>Sum+T (I +Iladd) 

NEXT I 

Tuin=Sum/6 

ELSE 

FOR I=0 TO 7 

IF (Ipos=0 AND I[*7) OR (Ipos=! AND I=6) THEN 3870 
Sum=Sum+T (1 +ladd) 

NEXT I 

Twin=Sum/7 

END IF 

PRINT USING “12X,""Avq@ Local Heated Wall Temp at Inlet (Tuin) no) 0) - 
( 


Cy ee huen 

! CALCULATE AVERAGE LOCAL HEATED WALL TEMPERATURE NEAR MIDDLE. 
Sum=0, 

! CORRECTION FOR BROKEN T/C CHANNEL 16 CURVED HEATER ONLY. AND T/C 


3961 !CHANNEL 11 STRAIGHT HEATER ONLY. 


S700 
Sis 513. 
3990 
4000 
4010 
4020 
4030 
4031 
4040 
4050 
4060 
4070 
4080 
4081 
4090 


IF Itype=! THEN 

FOR I=93 TO 16 

IF (Ipos=0 AND I=t0) OR (Ipos=0 AND I=!6) OR (Ipos2=! AND I=9) THEN 4020 
Sum=Sum+tT (1 +ladd) 

NEXT [ 

IF Ipos=0 THEN Tumid=Sum/6.0 

[F Ipos= 17 HE Tumra=Sum/ 7.0 

ELSE 

FOR [=3 10) 76 

IF (Ipos=0 AND I=ft0) OR (Ipos=! AND I=93) OR (Ipos=0 AND I=11) THEN 4080 
Sum=Sum+T(I+Iadd) 

NEXT I 

IF Ipos20 THEN Tumid=Sum/6.0 

IF Ipos=1 THEN Tuwmid=Sum/7.0 
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4100 
4110 
DD. ?" 
4120! 
4130 
4140! 
4150 
4160! 
4170 
4180! 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
Db 
4320 
4330 
4340! 
Ago 0 
4360! 
4370 
4380 
4390 
4400 
4410 
4420 
4430 
4440! 
4450 
4460! 
4470 
4480 
4490 
4500 
4510! 
4520 
4530! 
4540 
4550 
4560 
4570 
4530 
4590 
4600 
4510! 
4620 
4630! 
4640 
4650 
4660 
4670 
4680 


END IF 
PRINT USING °12X,""Avq Local Heated Wall Temp Near Middle (Twmid) = "",3D. 
bey "+ Tumid 


' CALCULATE AVERAGE LOCAL HEATED WALL TEMPERATURE AT QUTLET. 
Sum=0. 
! CORRECTION FOR BROKEN T/C CHANNEL 26 CURVED HEATER ONLY. 


IF Itype=! AND Ipos=0 THEN 
FOR I=23 TO 29 

IF [226 THEN 4230 
Sum=Sumt+T (1) 

NEXT I 

Twout =Sum/6 

ELSE 

FOR [#23 TO 29 
Sum=Sum+T(1+Iadd) 

NEXT I 

Twout2Sum//7 

END IF 

PRINT USING “12X,""“Avq Local Heated Wall Temp at Outlet (Twout) i opaaeereo) Iie 
CC) "> Iwout 
PRINT 

PRINT 


EHIER ORIFICE PRESSURE DATA. 


BEEP 

INPUT “ENTER PATM (CinHa), DPM(i1nH20), PIM(CinH20), REI". Patm,Dom,Pim.Rel 
PRINT USING “9X,""“THE FOLLOWING ORIFICE DATA WERE ENTERED:""" 

PRINT 

PRINT USING "12X,""Patm(inHg) DPM(inH20) PIM(1nH20) REiges 
PRINT USING "14X,2(2D.DD,10X) ,(2D.DD.7X) ,SD.D" ;Patm,Dpm,Pim,Rel 

PRINT 


! CONVERT PRESSURE READINGS TO SI UNITS AND P1M TO ABSOLUTE. 


Patm=Patm*3376.8 
Pdel=Dpm*248.34 
P1=Patm-(P1m#*248.84) 
Torf='(33) . 


! PHYSICAL PROPERTIES AND CONSTANTS. 


R=286.98/7 '!FOR AIR. 

Rho=P1/(R* (Tort +273.15)) ! DENSITY OF AIR AT ORIF PLATE. 
Gamma=1.40 !FOR AIR. 

Bae, !1Cp GOOD FOR Tblk BETWEEN 12 AND 33 (DEG C). 
Ge= 1. 

Sigma=5.669E-8 

Ecu=.!12 !EMISSIVITY OF COPPER. 


! DIMENSTONS OF CHANNEL. 


De=.006350 
Wid=.2540 
Ri=.297 
Pwet=2*(De+tWid) 
Ac=Dc+Wid 


5D 


4690 
4700 
4710! 
4720 
4730! 
4740 
4750 
4760! 
4770 
4780! 
4790 
4800 
4810 
4820 
4830! 
4840 
4850! 
4860 
4870 
4880 
4890 
4900 
4910 
4920 
4930 
4940! 
4950 
4960! 
4370 
4980! 
4990 
5000! 
Suu 


Dhd=4*Ac/Pwet 
Apl=.0774192 


! PROPERTIES QF INSULATION. 


Mine=.O 127 
Kins=3.8E-2 


! ENTER ORIFICE CONFIGURATION. 


BEEP 

INPUT “SELECT DIAMETER OF ORIFICE (0=.5334, 121.0755 (inehes))) sone 
IF Size=@ THEN Dorf=.013548 

IF Size=! THEN Dort=.027318 


! DIMENSIONS OF PIPE AND ORTR IGE aie 


A=(PI*Dorf°2)/4 
Dpipe=.051772 
Apipe=PI*Dpipe “2/4 
Beta=Dorf/Dpipe 


PRINT USING “12X,""Dorf(m) A(m*2) Beta”“"* 
PRINT USING “12X,(Z.6D.6X),(Z.3DE,6x), 2-4) = or feaeee es 
PRINT 
PRINT 


! CORRELATION FOR EXPANSION FACTOR BASED ON 1D AND 1/2D TAPS. 
Y=1-(€.41+.35*Beta’ 4)*(Pdel/(Gamma+P1)) 

! CALCULATION OF FLOW COEFFICIENT (K) IN STI UNITS. 

B= .0002+2.794E-S/Dpipet(.0038+1.016E-S/Dpipe) *( Beta 2+(16.5+1.968504E+2*Dp 


ipe)*Beta’ 16) 


S020! 
5030 
S040! 
5050 
S060 
+Beta’ 
S070 
S080 
S090! 
5100 
S110! 
S20) 
S130! 
S140 
SiS 
5160 
5170 
S180 
5190 
S200 
S210 
52c0 
S230 
S240 
6250 


CP aaa: 


! K1 AND K2 ARE DUMMY VARIABLES TQ CALCULATE Ko. 


K1=.6014-5.3974064E-3/Dpipe’* .2S 
K2=(.376+2.8971138E-2/Dpipe’.25)*(1.6129E-7/(Dpipe’2*Beta’2+6.35E-S*Dpipe) 
4+1.5*Beta’ 16) 

Ko=K1+K2 

K=Kot+1000*B/Re1l°.5S 


! Mu CORRELATION GOOD FOR Torf BETWEEN 17 AND 41 DEGREES CELSIUS. 
Mu=4,6971492Z9E-S*lorf*l. 71347 22b os 
! CALCULATE Mdot AND REpipe AND COMPARE [0 PREDTGtaoR ar 


Mdot=K+*A+¢Y¥+(2*Gc+*Rho«Pdel)*.S 

Repipe=(Mdot*Dpipe)/(Mu*Apipe) 

Dit f=(ABS(Rel-Repipe) /Rel)+*100 

Rel=Repipe 

IF Dit f>.001 THEN S080 

Coy USING "9X,"°"THE FOLLOWING DATA WERE CAUCUEATED = - 
N 


PRINT USING *12X.""Orifice Expansion Factor er) a OO OC Dineen 
PRINT USING "12X,*"Oritice Flow Coett teremcmue a 8h Os VOUS 
PRINT USING “12X.""“Density Based on Torf (Rha) = "tall (kq/m- 


Rho 


56 


5260 


woth 
9 


B20! 
5280 
ae! 
5300 
5310 
S320 
Sos! 
5340 
Bo50! 
S360 
5370 
5380 
Bo I0 
5400 
5410! 
5420 
5430! 
S440 
5450! 
5460 
5470! 
5480 
5490 
5500 
5510 
$520 
5530 
5540 
S500 
5560 
55/0)! 
5580 
Boo0! 
5600 
S610! 
5620 
6500! 
5640 
5650 
5660! 
5670 
5630! 
5690 
5700 
a7 10 
or 20 
a7 30 
S740 
= (sit 
57650 
a7 70 
S780! 
Sel 
5800! 
5810 
5820 


> Mu 
53830 


PRINT USING “12X,"“Viscosity Based on Torf (Mu) Ae SUE ee (OG) M.S) 


! Mu CORRELATION GOOD FOR Tblk BETWEEN 17 AND 41 DEGREES CELSIUS. 


Mu=4.6971429E-8*Tblk+1.7194/722E-5 
Red=(Mdot*Dc)/(Mu*Ac) 
Rehd=(Mdot*Dhd) /(Mu*Ac) 


! CALCULATION OF POWER INTO HEATER PLATE. 


ITadd=0 

IF Ipos=t THEN Iadd=39 

IF Ipos=0 THEN Rpr=2.01729 

IF Ipos=! THEN Rpr=2.00839 

Qp=(Emf (38+Iadd) *Emf (39+Iadd))/Rpr 


! CALCULATION OF HEAT INPUT TO AIR. 
Qair=Ndot*Cp*(Tout-Tin) 
! CALCULATE LOCAL HEAT TRANSFER COEF. AND NUSSELT NUMBER. 


Kair=/7./25/7143E-S#Tcint+.024165836 
Hin=Qair/(Apl*(Tuwin-Tcin) ) 
Nuin=(HintDhd)/Kair 
Katr=/.7257 143E-5* lemid+.024165836 
Hmid=Qair/(Apl+*(Tumid-Tcmid) ) 
Numid=(Hmid*Dhd) /Kair 

Kair=/.7257 143—E-S<«Tcout+.024165836 
Hout=Qair/(Apl*( Twout-Tcout) ) 
Nuout=(Hout*Dhd) /Kair 


! Kair CORRELATION GOOD FOR Tblik BETWEEN 17 AND 41 DEGREES CELSIUS. 
Kair#/7.725/7143E-S*Tb1k+.024165836 
! CALCULATE AVERAGE HEAT TRANSFER COEF. AND NUSSELT NUMBER. 


Havag=Qair/(Apl*Tdel) 
Nuavg=(Havg*tDhd) /Kair 


Me cobCULH iE wien) IU0SoES. 


Qlol=(Kins*Apl*(Tinsol-Tinso2))/Xins 
Qlo2=(Kins*tApl*(Tinso2-Tinso3))/Xins 
Qlo=(Qlol+Qlo2)/2 
Qlil=(KinstApl#(Tinsil-Tinsi2))/Xins 
Qli2=(Kins*Ap] +(Tinsi2-Tinsi3))/Xins 
Q11=(Q01114+0112)/2 
Rr=2.0+((1l-Ecu)/(Ap|+Ecu))+(1.0/Ap1) 
Qr=(Sigma*(ABS(( Two+273.15) °4-(Twit273.15)°4)))/Re 
Qdel =Qp-Qair-Qlol-Qr-01:1 


{! CALCULATE DEAN NUMBER. 


De=.5*Red*#((Dc/2)/R1)*.5 
PRINT USING “12X,""“Viscosity Based on Tblk (Mu) Ser ee ODEs = Ck avin. S.) 


PIU SING= 12k. iherm Cond Based on Iblk (Kair) = ““,Z.3DE.°" (W/m.K)" 


3) 


© Kaiten 


53840 
S850 
5860 
53870 
$380 
5890 
5900 
S210 
Soo! 
$930 
S340! 
eed) 
5960 
5370 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


USING 
US ING 
USING 
USING 
USING 
USING 
USING 


12k. atm 
"Leky: del 
alta hee 


Hens sandot 


"12k." "Reptpe 


"12X," "Red 
"12X,""Rehd 


! ONLY PRINT DEAN NUMBER 


IF Itype=1 THEN PRINT USING “12X,""De 

PRINT 

PRINT USING “12x,”""Qp EN oO Da 
atte)" 0p ,Gain 

PRINETCUSING= “12%. ""Glol =. “* 2220220. 


5380 


(Watts) >>" <Glol-Gin 


$990 


PRINT USING 


ZK y, me oZ 


(Watts)? ™“SOlez 0 


6000 


(Watts) "<O0le0ly 


6010 


(Natts) "Or cael 


6020 
6030 


PRINT USING 
PRINT USING 


PRINT 


PRINT USING 


in. Nuin 


6040 


PRINT USING 


mid.Numid 


6050 


PRINT USING 


out,Nuout 


6060 


PRINT USING 


avg. Nuava 


6070 
6080 
6090 
6100 
5110 
6120 
6130 
6140 
6150 
6160 
6170 
61380 
6190 
6200 
Seat 
6220 
6230 
5240 
6250 
6260 
6270 
0280 
6290 
6300 
6310 
6320! 
6330 


— 


PRINT 
PRINT 
BEEP 

INPUT 
IF Go_on<>0 
PRINT USING 


lek, . role 


Sik Ue 


Sick in 
"V2, -ohimed 
ek. : howe 


“12X,""Havg 


THEN 
ssa ae a 


IF Im=0 THEN WAIT 300 


GOTO 1700 
END iP 


PRINT USING 
PRISE Roe 


al da 


ASSIGN wFile TO * 


END 


o> ec De oa 
="! 2D. 3D ae 
o Ge Dees, 


SU) UD iw 
. 730200 
2 OUD 
m4 30. DD, ** 


“PRINT USING “3X7 “END OF RUNG es 


|. CONVERTS EHF 10 DEGREES GER sIGS 


SUB TvysviV,T) 
COM. /Go7 D7 aati.) 


Sum=0 


FOR I*0 70 7? 


Sum=Sum+D( 1) #V "I 


NEXT I 
T=Sum 
SUBEND 


! CALIBRATES 1/C READIUGS. 
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6D .DDs ee 

oY AD. DD 
“* 6D Obes 

ee eer Die 
"".6D.D";:Repipe 
"* 60.D" Red 
"" 6D.D"3;Rehd 


(Hatts)°"* «7X Carne 


(N/m°2)""":Patm 


(N/m 2G so ded 


(N/m 2) en 


(ka/s)"""sMdot 


WHEN HEATING CURVED SECTION. 
- "",4D.D";De 


(Watts) 7°". 7) 40 aie 
(Watts) "" 7%. > 20 ieee 


(Watts) ” 72 20) tie 


(Watts) ° 7.7%. Ode 


(W/m 
(W/m 
(W/m 
(W/m 


"2C)"" 0X 
"2C)* "7X, “Numed = 
"20"? (7h. NuoU tees 
"2C)"",7X,""Nuavg = 


“WILL THERE BE ANOTHER RUN? (1t=Y,0=N)*,Go_on 


“* 3D.30 ae 


ere 


ww eDe 
pn Ne 
wry 4)b): 


hes) 218) 
ee Ue 
> our 
“ee 0) 


(W 
305 
30 
3D 
SD ya 


DD,"<H 
DD, "3H 
DD, "3A 
DD,“+H 


6340! 
6350 
9560 
6370 
6380 
6390 
6400 
6410 
642 


Vereen iecors (i. ! ) 

COlMey Coz + D7). Aal7Sec) 
Tte=Aal lO) 

FOR J=1 10 2 
To=TotAal(I,J)*T*J 

NEXT J 

RETURN Tec 

FNEND 


Sy, 


APPENDIX B 
SAMPLE CALCULATION 


Figure (B.1) shows the components of the heat balance in the test channel. A 
sample calculation follows that details the methods used in the data reduction program 
to estimate the value of the energy balance components and then calculate the 
hydraulic Reynolds number and average Nusselt number for each data set. A copy of 
an actual record of one set of experimental data for the concave test section precedes 
the sample calculations together with a list of the parameters necessary to make the 


sample calculation. 


HEATED COPPER WALL 


Qlo 





Figure B.1 Test Section Energy Balance. 
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1. SAMPLE CALCULATION DATA 


hen th. 
THE FOLLOWING DATA WERE RECORDED: 
91.7554) 


Heater Voltage (Vh) 
Resistor Voltage (Vr) 


Pree toon Reslotanmnor? 
GITeRePLale: 
ie eetlimner: 1 it 
heme 0G) : fee os ech ylly 
le Cahlamie pr: 7 o 
henewee! rps) esis Ya a3 
T/C onber: 13 14 
Tememcc i; TNS Sy 
T/C Number: bs 20 
lene UG: 1 Ses 84.26 
T/C Number: oS 26 
Teme (0): (hi eye Slee 
OUTER INSULATION: 
/C Number: ei Clinso |) 
Tempe te): SZ eae 
ORG IGE Mien st lortye=. 30.25 CO) 
INLET TEMPERATURE: 
T/C Number: 2S 36 
Lemon (eG): e703 22 e 
INNER PLATE: 
T/C Number: 4} 42 
Temp (0); 23.89 Oe es! 
76 Number: 47 48 
Teme. (GC): Coace Cong 
T/C Number: 53 C4 
femoece): od | Sa 3 
T/C Number: £0 
temp (CC): 23.80 24.49 
17C Number: 65 66 
Temp (C): 25.56 Chae ae ey 
INNER INSULATION: 
C Number: fie tl insa |? 
Temp (C): 24.401 
OUTEET TEMPERATURE: 
T/C Number: 74 ios: 
Temp (C): SoG Sats 


Day. and Time: 


Hh nw tt 


6] 


7.704 (V) 
Ces AOinms) 


eevee 4.54 


ci 4 
76.00 68.01 
5 10 
(eles. (dob 
Ks 16 
Vs 0 Boao? 
a oe 
Sci 69.58 

O37 fi 
47.70 73225 
2 (linse2) 

56.88 

of ais! 

e208 eo 07 
3 44 
23 240 Poo 

49 S0 
23.56 Powe 

oS S56 
2238 Congo 

61 b2 
Cae, Cuezor 

67 68 
eComlic 2f04 

Chencic) 
Pho 

76 ai 

30.65 30.69 


(Jlinso3) 


40.01 


THE FOLLOWING TEMPERATURES WERE CALCULATED: 


30.32 30 
73,62 °NGy 
eval, WC, 


Average Local Channel lemp at Out lerrt leone) 
Ava Local Heated Wall Temp at Inlet (Twin) 
Ava Local Heated Wall Temp Near Middle (Twmid) 


Average Quter Wal! Temperature (Two) = 76.46 (C) 

Average Inner Wail Temperature (Tw) = 24,15 (0) 

Average Outlet Temperature (Jout) = 3060 meg 

Average Inlet Temperature (Tin) = 22.08 

Channel Inlet and Outlet Temp Diff (Tdiff) 2 3.61 (Gd 

Average Bulk Temperature (Tblk) = 26.35 sce 

Mean Temperature Difference (Tdel) > SO -tsee 

Outer and Inner Hall Temp Difterence (Tudiff) ~=* Pate) (C0. 

Average Local Channel Temp at Inlet (Tein) = 22.44 (C) 

Average Local Channel! Temp at Middle (Temid) = 25,3076) 

Avqg Local Heated Wall Temp at Outlet (Twout) = (9273 50e) 
THE FOLLOWING ORIFICE DATA WERE ENTERED: 

P3tm( inHg) DPM(1nH20) P1M(1nH20) RE 1 

29466 17.50 Po. 30 1500020 
Dorf (m) AGm 2} Beta 
O07 onre SoG le v4 0. Save 


THE FOLLOWING DATA WERE CALCULATED: 


Ortfice Exeansivon Factor a. ~ 0.959 

Orifice Flow Coettierent (Kk) =~") OS 

Density Based on Tort (Rho) oS (ka/m~> 3) 
Viscosity Based on lor remains = 18.616E-06 (kaq/m.s) 
Viscosity Based on JIbIe (My) = 13.434E-06 (ka/m.s) 
Therm Cond Based on Tblk (Kair) = 2.620E-02 (W/m.K) 


Patm = 100848.13 (N/m*2) 

Pdel = 4579256 Clyne) 

Pal = 96642.74 (N/m*2) 

dot = 0.0360 (ka/s) 

Repipe = 4/7611.8 

Red = 76072 

Rehd = TS0S4 

De = Gis Fares 

Qp = 350.417) CWatts) Gair = Stl. dol Gates) 
Qlol = 5.308 CWattse ie = 063 (Watts) 
NQlo2 = 5, 008 (ieaets) By ys » LOOM Cha tie. 
Whos <= 4,308 (Ha crs3 lio = mbes Clee a7. 
Or = [296 SWati=) del = 30.453 GNates) 
Hitee=" Ged 4> Cy meetu Mutim ~ = 637s 

Hmid = 64.38 (W/m 2C) Humid = 40.03 

Hout = 81.48 (W/m’> 2c) Nuout = 38.08 

Havg = 80.47 (W/m'>2C) Nitaval =" ous 


END OF RUN 


<< 
Ss 


< 
5 


ac) as 
ae} 
Lees | 


atm 


te 
22 


meee 


rm - TQ 


Q 


TW 


WO-Wil 


pl 


c 


> > 


ae 


SE eS 


pipe 


orf 


0 


= 7.704 V ie 


le o5 V IP eee), 11) 


welt 
mn 


2.01729 Q 
29.86 in. Hg 
17.60 in. H,O 
16.90 in. H,O 
POS C 

30.68 °C 

76.46 °C 

Decne @ 

Deron 

pe a 

Oe © 

82.38 °C 

56.88 °C 

40.01 °C 

30.25) Te 

0262 W:m-°K 
1007 J Kg-°K 
4.6971 X 10° x T.. + 1.71947* 105 Kg/m-s 
.0380 W/m-°C 
0127 m 

0.12 

5.669 X 108 W/m?-°K4 
Soy 

1.40 

NOUK Gems 
286.98 Nm, kg-°K 
1.0 

07742 m? 

0016 m? 

00211 m? 

00635 m 

O518 m 

027318 m 
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2, TEMPERATURE CALCULATIONS 


a. Average Temperature of the Outer (Inner) Wall 


y Outer (Inner) Plate Temperatures 
T = 76.46°C (24.16°C) 


wol Twi) — Number of Outer (Inner) Plate Thermocouple 


b. Average Channel Outlet (Inlet) Temperature 


¥’ Outlet (Inlet) Temperatures 
(T..) = — = 3 30.68°C (22.08°C) 
our in” Number of Outlet (Inlet) Thermocouples 


c. Channel Inlet to Outlet Temperature Difference 


Var = lou lian > OU? 6 ee Cee 

d. Bulk Temperature of the Flow Through the Channel 

ences | 2205 © Pees 2 
1. -  S — ore 
os 2.0 2.0 

e. Average Temperature Difference of the Heated Wall and Fluid Bulk Temperature 

AT = T... - Ty, = 76.46°C - 26.38°C = 50.08°C 
WO blk 

f. Temperature Difference Between the Outer and Inner Walls 

Vers ee ae a OO GG ee 


3: POWER CALCULATIONS 


a. Power supplied by the Heated Plate 


om Mn Mor. 21755 477 
P Ra, 2.017292 


are |} 


50.41 Watts 
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b. Power Loss Through the Outer Wall Insulation 
Qi ~ (1/2)Quo) ic Qio2) 


= (1/2){( Er ago ee yk. 
AX... Ins A, 


a (—inso2 : cos ‘eek Me 
ins’ pl 
““ins 


(0.038 W'm°C) 


Weegee 182938 56.88(°C)) + (56.88 - 40.01(°C))] = 


LO sFioal 


c. Power Loss Through the Inner Wall Insulation 


Q3 = C2KQ + Qi) 


= (}/? Wage - Ti nsi2 Le 
<< ( j 2){( ee ) ins An 


p tins“ Tins3_ ee 
( AX. ) ins a 
Ins 
5 (0.038 W'm’C) 
(1/2)(.07742m i) aa a aoe C)) ce 24.13 23.90( ©) = 
LO Fane’ 


d. Heat Loss by Radiation 


l. Radiation Resistance 








- € l 
R. = 32 sc sre re 
Ay : eet An . ae Wi 
l 2 l 7 
An fou 0.07742m.. iw 12 
2. Heat Radtated 

4 4 

Q = os ak i 
r iS 

(5.669 x 10 §)[(76.46+ 273) - (24.16 + 273)4] 


= 1.996 W 
202.36 


4.908 W 


O81 W 


4. MASS FLOW RATE CALCULATIONS 


a. Pressure Conversions 
Pim = 29-865 in Hg x 3376.8 = 100848.13 Nm? 
AP = 17.60 in H,O X 248.84 = 4379.58 N’m? 


P, = (16.90 in H,O x 248.84) - 100848.13 = 96642.74 N'm? 


b. Density of Air 
(96642.74) 


c. Expansion Factor 
AP 
Y= 1-O8l + 03507) = 
a 
4379.58 


= | - (0.41 + 0:35(0.5277)) —$_$_—$—— roe 
(1.40)(96642.74) 


d. Mass Flow Rate 
hh = YKA /2gp,, AP = (.9859)K(.0005861m1’) 
2(1 Kg-m, Ns”)(1.1099 Kg m*)(4379.58 N’m*) = .05669K Kg's 


Iterate to find m: 





Assume a hydraulic Reynolds number = 15000 


: Deine cy pe =  O5ISX00I6 
“pine ~ (Da) Sha ~ “Oizzgx 00211) > > “hd 
Pip 


+ Rei... ~ 15000(3.1961) = 47942 


Interpolate in the appropriate table in Reference 15 for K, the flow coefficient: 
kK = .6329 


(Note that the computer program used to reduce the data used an approximating 
polynomial to find the correct value of the flow coefficient, K. 


> m = (.05669)(.6329) = .0359 kg s 


Solve for the improved value of er ee 


ih Die” _1(00359)(0; 1c) 
Soe. I (.00211)(18.616 X 10°) | nee 
; pipe! air ee 


Iterate until successive values of m differ by less than .001: m = .0360 
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5. REYNOLDS NUMBER CALCULATIONS 


a. Re Based on Pipe Diameter 


mD.. .0360).0518 
Re eepipe: = Seah = 47475 
pipe A op (.00211)( 18.616 X 10 9) 


pipe’ air 


b. Re Based on Channel Height 
ee) te as (.0360)(.00635) 


hd 
A 


air 


6. HEAT TRANSFERRED TO THE WORKING FLUID 


Qair = MC, (Thy, - Tig) = (0360 Kg, 1007 J/Kg-°C)(30.86 - 22.08)°C = 312 W 


7. AVERAGE HEAT TRANSFER COEFFICIENT CALCULATION 


AAT (.07742m*)(50.08°C) 


8. AVERAGE NUSSELT NUMBER CALCULATION 

; (80.5 W/m*-°C)(.01229m*) 

Nu, = —2*- = ——______——_ = 378 
kK. 10262) Worn © 


9. DEAN NUMBER CALCULATION 


2 7751 /.00635/2 


297 
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APPENDIX C 
EXPERIMENTAL UNCERTAINTY 


The experimental uncertainty in calculated parameters was determined in 
accordance with the method described by Kline and McClintock [Ref. 23]. 


The qualitative results shown here correspond to the data set and calculation of 
Appendix B. There was a sigmificant decrease in the experimental uncertainty of the 
average Nusselt number for Reynolds numbers above 17,000 in this experimentation 
compared to that of Galvo [Ref. I4: pp. 104-105] because of the changes made to the 


test channel that promoted mixing of the working fluid at the exit of the channel. 


The equations used to determine the experimental uncertainty of calculated 


parameters are shown below. 


1. DENSITY OF THE WORKING FLUID 


=e —— 4 ~_ Ss -~ ao cae) abet oo ‘ 
; [( 2 ) ( 5 ) ( 7 me 


l orf 


2. MASS FLOW RATE OF THE WORKING FLUID 


mm Yo. Ky Ap, 
rh ae K ae 
$0. ANS 
Plat (eee ro ey) Oe 

pag AL 


3: HEAT TRANSFERRED TO THE WORKING FLUID 








| 1 oe 
(Ore 2 (By (AY 
Ge a 
¢ ee is 
Gree an) 


6S 





4. AVERAGE HEAT TRANSFER COEFFICIENT 
oh = | OQ, ? + ( OA )? + ( oa es 
h Ay Nel, 
air p 
5. AVERAGE NUSSELT NUMBER 
6Nupg _ , Sh, 8Dng 2 4. OKaie \2 10.5 
= [( r J uaa ) ( z ] 
aad hd air 
6. HYDRAULIC REYNOLDS NUMBER 
dRe = [( ils )? ue ( ODig ? a ( OW, )? a ees \- be 
R mD hd Mair A 


Several parameters used in the calculations described herein had uncertainties 
that were constant over the entire range of testing. These quantities and their relative 


uncertainties follow. 


Quantity (B) Uncertainty (+06B,B) 
A 0020 
A. 0408 
Ay 0006 
© .0020 
D. 0400 
Di 0391 
K 0030 
Kee 0038 
a 0010 
R 0003 
Wi 0080 
bd 0020 
L.. 0054 
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The approximating polynomials used to resolve thermocouple voltages into 
temperatures matched actual temperatures to within less than+0.1°C over a range of 
temperatures from 20°C to 85°C; however, in some cases this temperature range was 
exceeded. A conservative estimate of the uncertainty of temperatures measured by 
thermocouples was set at £0.1°C, the accuracy of the thermocouples specified by the 


manufacturer of the thermocouple wire. The following uncertainties were used in 


calculations: 
Quantity (T) Uncertainty (= 
ie OG 
ie Ole 
J 0.1 °C 
orf 
or: Orlane 
ler Os. Ce 
a ee 0.236 
AT OS 


For the data taken at a hydraulic Revnolds number of 15000 in the concave test 


section the following additional uncertainties applied: 


Quantity (P) Uncertainty (7) 
i 0.05 
AP 0.10 


The results of the uncertainty calculations at a Reynolds number of 15000 in the 


concave test section were: 


Quantity (B) Uncertainty (+6B/B) 
p 00445 
m 00544 
Or 02397 
h 02431 
a 04620 
Rew 05703 


The major contribution to experimental uncertainty in both the average Nusselt 
number and hydraulic Reynolds number was the uncertainty in cross sectional area of 
the channel. The hydraulic diameter of the channel was a function of channel cross 
sectional area, resulting in a relatively large uncertainty in the hydraulic diameter. 
These two parameters accounted for approximately fourty percent of the uncertainty in 
the average Nusselt number calculation, and about eighty percent of the uncertainty in 


the calculated hydraulic Reynolds number. 


TI 


o>) 


Can 
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